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Estimating the fidelity of state preparation in multi-qubit systems is generally a time-consuming task. Never- 
theless, this complexity can be reduced if the desired state can be characterized by certain symmetries measur- 
able with the corresponding experimental setup. In this paper we give simple expressions to estimate the fidelity 
of multi-qubit state preparation for rotational-invariant, stabilizer, and generalized coherent states. We specifi- 
cally discuss the cat, W-type, and generalized coherent states, and obtain efficiently measurable lower bounds 
for the fidelity. We use these techniques to estimate the fidelity of a quantum simulation of an Ising-like inter- 
acting model using two trapped ions. These results are directly applicable to experiments using fidelity-based 
entanglement witnesses, such as quantum simulations and quantum computation. 
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I. INTRODUCTION 

Highly entangled states provide required resources for 
quantum information processing (QIP), a developing field ad- 
vancing both the fundamental understanding of quantum sys- 
tems and novel technologies. Entangled states are used to en- 
code qubits for fault-tolerant quantum computation jT| and 
for more efficient quantum state readout [Zj. Entangled states 
are used for quantum communication over long distances and 
teleportation protocols ||3|. Finally, highly entangled states are 
central to many-body quantum simulations, whose power lies 
in their ability to coherently manipulate such states for later 
analysis j^QSQl- Entangled-state preparation in any QIP 
system, and its verification, is thus of paramount importance. 

One successful architecture for QIP is the trapped-ion sys- 
tem, in which qubits are encoded in the internal electronic 
states of ions, and laser fields can control the collective inter- 
nal and external states of the ions. Recently |(^ j^il, multi- 
qubit entanglement has been experimentally demonstrated 
in these devices. In Ref. J^l, quantum state tomography 
(QST) I id fTlll was employed to verify that W-type states for 
up to = 8 ions (qubits) were produced. Since the dimen- 
sion of the Hilbert space H associated with a quantum system 
increases exponentially with the system size (as does the di- 
mension of the density matrix), performing full QST is, in 
general, extremely inefficient for large systems. For exam- 
ple, realizing QST on an ion-trap device requires on the order 
of 0{3^) measurements, where N is the number of qubits 
involved that are measured in the x, y, and z-bases. In 0] 
the full QST process for iV = 8 ions required 656, 100 mea- 
surements over ten hours. This extremely large data set re- 
duced errors due to quantum projection noise 1 12], until other 
sources of error (such as imperfect optical pumping, ion ad- 
dressing errors, non-resonant excitations and optical decoher- 
ence) dominate. Such examples illustrate a potential road- 
block to practical implementation of large-scale QIP: it is im- 
possible to exploit the speedups associated with QIP if an 
exponentially-large amount of processing must be performed 



to verify the creation of the desired states. 

It is important then to investigate efficient methods to es- 
timate the reliability of experimental quantum state prepara- 
tion. Here we point out that many useful entangled states have 
certain symmetries which allow fidelity determination without 
full QST. For these states, an efficient number (polynomial in 
N) of measurements is sufficient to obtain lower bounds for 
the fidelity. A similar technique has been used to determine 
a lower bound on the fidelity of several -particle cat states UJ; 
we describe and generalize such methods. To see this, we 
use the quantum fidelity as a measure of the distance between 
quantum states Specifically, the quantum fidelity be- 
tween the actual state prepared in the laboratory pi, which is 
in general mixed (i.e., Tr{pf) < 1), and the desired pure state 
lip) to be prepared is defined by 



(1) 



Equation Q can be evaluated by measuring the expectation 
value of the density operator p.^ = over the state pi. 

For example, if is a product state, then p^, has only one 
non-zero matrix element (in the right basis) that is along its 
diagonal. The fidelity !F{pi,p^,) can be simply obtained by 
repeatedly preparing pi and then measuring the population of 
the state \tp). 

More generally, the density matrix of an iV-qubit system 
is a linear combination of operators belonging to the u(2^) 
algebra: 
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where the subscripts aj = 0,1,2,3 correspond to the Pauli 
operators 1, ax, cry, and az, respectively. (The symbol ® rep- 
resents the matrix tensor product.) These operators are given 
by 
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In particular, cr^ . = 1^ ® • 
with the Pauli matrix ctq being located at the jth position 



2 



in the decomposition. From now on, we remove the symbol 
(g) from the products of PauH operators. We also adopt the 
convention {A)p^ = Tr[p(A]. 

are given by 



The real coefficients Cai,a2, 



= 2 



-N 



= 2-^Tr[p(ai^...a^„)] (Otherwise) 



Then, 



(4) 



(5) 



in terms of these operators only. Thus, the fidelity of having 
prepared [Eq. Q] can be estimated by measuring observ- 
ables, over the actual prepared state pi, that solely involve cor- 
relations between the O^'s. In other words, measurements in 
bases not related to the symmetry operators are not required 
because they do not provide any information when evaluat- 
ing the fidelity of state preparation. The purpose of this sec- 
tion is then to give lower bounds for estimating the fidelity 
of state preparation for three classes of highly-symmetric A^- 
qubit quantum states, and show that these can be efficiently 
obtained. 



and full QST is generally needed to estimate the coefficients 
Cai, ,Qjv required to evaluate Eq. (|5}- However, if the state 
\ip) can be uniquely characterized by certain symmetries, 
some of the coefficients Cq"^,....q„ will vanish and the cor- 



responding 



need not be measured. Full QST over 



pi is then no longer required, and the complexity of evaluat- 
ing Eq. (|5} or of setting a lower bound on J-'^{pi, p^) can be 
greatly reduced. 

A straightforward example of using symmetry to sim- 
plify fidelity estimation can be seen in previous work 
with iV-qubit cat states \Q\\Z)m = ^(|0i02 • ■ • Oat) + 
II1I2 • • • Iat)) in trapped ion systems fl3l H . The 
I GHZ) TV state is uniquely defined by the symmetry opera- 
tors {cricrl ■ ■ ■ cr^, cr^c^, cr^cr^, • • • , a^'^cr^}, that leave the 
state unchanged after their action. As we will show, the fi- 
delity of having prepared |GHZ)jv can be estimated by mea- 
suring the expectation values of the symmetry operators. In an 
ion-trap setup, for example, repeated simultaneous measure- 
ments of the projections of all of the the ion spins along the 
X axis, and of all of the ion spins along the z axis, gives the 
fidelity of having prepared the | GHZ) at state 1 15]. 

In Sec. in] we expand this idea to study certain cases in 
which the desired state can be characterized by different types 
of symmetries. First, we focus on the class of rotational- 
invariant states (i.e., eigenstates of the total angular mo- 
mentum operator) since some interesting entangled states for 
quantum information tasks are in this class L16J. Second, we 
study the family of stabilizer states (SSs) which provide the 
foundation of the stabilizer formaUsm used in different quan- 
tum error-correcting procedures iITtIi . Third, we study the case 
of generalized coherent states (GCSs) which provide a natu- 
ral framework to stud^ertain quantum simulations of many- 
body problems Q B- In Sec. |in| we apply the obtained 
results to estimate (numerically) the fidelity of evolving the 
internal states of two trapped ions with an Ising-like Hamil- 
tonian, using the methods described in Ref. Finally, in 
Sec. lIVI we discuss the estimation of the fidelity of state prepa- 
ration due to the statistics from a finite number of experiments, 
and in Sec.|V]we present the conclusions. 



n. QUANTUM FIDELITY AND HIGHLY SYMMETRIC 
STATES 



The density operator of a pure state | ?/') , uniquely character- 
ized by its symmetry operators {Oi, ■ ■ ■ ,0l}, can be written 



A. Rotational-invariant states 

For a system of N qubits, the rotational-invariant pure 
states are completely specified by the equations 



' jz IJi jz 



(6) 
(7) 



where J'^ = + Jy + is the (squared) total angular mo- 
mentum operator, Jj = a}^ + + ■ ■ ■ + (7 = x, y, z), 
and (T^ is the corresponding Pauli operator acting on the jth 
qubit. The factor 2 in Eq. (|6j is because we are using Pauli 
operators instead of the actual spin- 1/2 operators. Then, the 
quantum numbers j and jz satisfy the following properties: 
j-max ^ j-max ^ |^^-| > 3, | Aj^ | > 2, and -N <j,<N 
(the symbol A indicates the difference between the corre- 
sponding eigenvalues). In particular, if —iV + 2 < jz < N — 2 
the state \j, j,) is entangled and for j = N, j^ = N — 2, then 
|A^,A^-2) = \Wn), v/ith 



\Wn) = ^[|li02---0jv) 
V N 



|0il2---0jv) + 
|0i02---ljv)]. (8) 



Although the \Wn) states are not maximally entangled for 
N > 2, they are particularly useful for processes such as tele- 
portation ll20ll . 

The density operator pjj^ of an iV-qubit rotational- 
invariant state with quantum numbers j = N and —N < j^ < 
N, in terms of the symmetry operators J and J^, is 



n 



-j<ji<j 



n 



0<j"'<7V 



(9) 



where tt^. = (J^ - j'^) and %.y = [J^ - j'{j' + 2)]. The 

symbol W denotes that the term has been excluded from 
the product. The normalization constant n is given by 



n_«,<,n 



b.-jOb-0- + 2)-.7'(/ + 2)]. 

(10) 

To evaluate the fidelity of Eq. Q, that is T{pu pj,j^) ~ 
[Tr(/5/pj j^)]^/^, it suffices to obtain the expectations of the 
correlations between the operators and appearing in 
Eq. (|9jl only. Although this procedure is still inefficient and 
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an exponentially large number (with respect to N) of observ- 
ables (i.e., products of Pauli operators) must be measured, it 
is more resource-efficient than performing full QST to obtain 

For example, if one is interested in preparing the Bell state 
I Bell) = \j = 2, J, = 0) = ^[|li02) + |0il2>] on an ion-trap 
device, the fidelity of faithful preparation could be obtained 
by performing measurements over three different bases only, 
corresponding to the expectations (cr^.(T^)p,, (o'yCry)p,, {crDp^, 
{al)p,, and {alal) p^,x&&^ecme\y. 

To obtain a lower bound on the fidelity of rotational- 
invariant state preparation, for j = N, we first define the oper- 
ators 5,/^. = -l{J,-j,)^andSj2^-^{J^~N{N + 2)). 
These satisfy 

with Gjvjv < -1 for ^ {j,jz) and e^j^ = 0. There- 

fore, for a general pure state = '^j' J'JJ' ^ Jz)' we 

obtain 

(12) 

where |cj.jj^ is the probability of projecting 10) onto the 
state \ j — N, jz) (i.e., the squared fidelity between the states). 
Since the actual prepared state pi is in general a convex com- 
bination of pure states, Eq. il2\ yields to 

^HpuPj.jJ > {Sj^+Sj2)p,+1. (13) 

This lower bound can be efficiently estimated by measuring 
only the observables J^, and a large number of times 
over the state pi. This corresponds to the measurement of 
3N^ — 2N expectations of different products of Pauli opera- 
tors; that is, polynomial in N. 

When j < N, the subspace with quantum numbers j, 
is degenerate and Eqs. (|6jl and Q do not specify the state 
uniquely. Then, Eq. (|9} becomes the projector onto the 
corresponding subspace. Nevertheless, the squared fidelity 
^"^{Pu Pj-j,) will still denote the probability of having cre- 
ated a pure or mixed quantum rotational-invariant state with 
quantum numbers j, j^. When j < N, the operator 5/2 must 
be redefined as 5 72 = — -^(J^ + 2))^, so the properties 
for the coefficients ej/_j/ in Eq. il 1> still hold. In this case, 
a lower bound to !F'^{pi, Pj.j, ) can be obtained by measuring 
iV'*/4 - iV^ + lliV2/4 - n(n even) expectations of different 
products of Pauli operators. 

B. Stabilizer states 

Another interesting family of states are the stabilizer 
states |[17il . which are defined by 

aiV^) = +l|V^) ;.se [1,5] . (14) 

The stabilizer operators Os e ii(2^) are products of Pauli 
operators fJlj and have ±1 as possible eigenvalues. (Note that 



Oi = l is the trivial stabilzer) An immediate consequence of 
Eq. ( I14> is that the operators Os commute with each other: 
[Os, Os'] = 0. Here, we focus on the case when the state |?/') 
is uniquely defined by Eq. (I14> : that is, the dimension of the 
stabilized space is one. The set Gs = {Oi, • • • , O5} forms 
the so called stabilizer group for Itjj). For practical purposes, 
we define Og in a compact way by its L linear independent 
generators lO: Gs = ■ • , satisfying 

g,\^) =+m; ie[l,L]. (15) 

Without loss of generality we can write = 

l5i = 1, ■ • • ,5l = !)■ 

The eigenstates of the stabilizer operators (associated with 
the stabilizer state) form a complete set of the 2^ dimensional 
Hilbert space H. Therefore, the density operator p^ can be 
written within this formalism as (1 = 1^ (g) • • • (g) 1^): 

Pi> = I51 = l,---,5L = = i,---,5L = i| 

1 ^ 

1=1 

and the fidelity [Eq. Q] can be estimated by measuring, over 
the actual state pi, the expectations of operators appearing in 
Eq. (d. 

A lower bound on the fidelity can be obtained in this case 
by defining the operator Sqs = U(Ylf=i9i) ~ {L ~ 2)11]. 
Then, 

'^Gslsi, • • • ,.9l) = egi,...,g^|gi, • • • ,5i>, (.9i = ±1), (17) 

with ei,..._i = 1 and Cgj^ ... < otherwise. Following the 
same procedure used for rotational-invariant states, we arrive 
to the inequality 

^^ipi,Pi,)>{SGs)pn (18) 

which can be efficiently estimated by measuring the expecta- 
tions {g,)p, Vi G [1,L]. 

As an example we consider the the Bell state |Bell) = 
-^[|0il2) ^ II1O2)]. For this state, the stabilizer group is 

defined by the generators Gs = {—ala^, —alal). Then, 
i = 2 and Sgs = ~ '^x'^x]- Another example 

is the set of maximally entangled 7V-qubit states |GHZ)Ar = 
-^[|0i02 ■ • - Oat) + II1I2 ■ • ■ Itv)]- For these states, the gen- 
erators of the corresponding stabilizer group are given by 
Gs = {cT^al • • • , alalajal • • • , - ), as pointed 
out in Seed For = 3, L = 3 and Sgs = W'^l^l^^l + 
^Wz+(^Wz-% 



C. Generalized coherent states 

The last class of states we consider are the generaUzed 
coherent states (GCSs) 12^1 . For a semi-simple, compact, 
7\ /-dimensional Lie algebra f) — {Qij Q2, • ■ ■ , Qm}, with 
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Qj = {QjV the iV-qubit operators acting on the 2^ dimen- 
sional Hilbert space H, the GCSs are defined via 

|GCS) = e'f'lhw). (19) 

Here, e''' denotes a unitary group operation (displacement) 
induced by (): e*'' = exp[i{J2j >^jQj)], G M. The 
state |hw) is the highest-weight state of f). To define it, one 
needs to assume a Cartan-Weyl (CW) decomposition t) = 
()n ® f)+ e [)" The set ()d {hi,-- - ,hr} is the 

Cartan subalgebra of f) (CSA) constructed from the largest set 
of commuting operators (observables) in [). The weight states 
\(j)i), which form a basis of states for H, are the eigenstates of 
i)D: 

hk\^,) = ul\(j),), k e [l,r], z e [0,2^ - 1] . (20) 

The sets ()+ = {e+ , • • • , e+ } and f)" = {e"^ , • • • , e" } are 
built from raising and lowering operators (e+ = e~t), and 
either map weight states into orthogonal weight states or an- 
nihilate them. (The subscripts Oj G W are the roots of f) and 
are considered to be positive.) Then, |hw) is defined by 

hk\hw) = Vk\hw), ke [l,r], (21) 
e+ |hw) = 0, (22) 

with Vk = u1 (i.e., we have assumed |hw) = \(j)o)). Note that 
M = r + 21. In many cases, |hw) = IO1O2 • • • Ojv), where \0i) 
represents an eigenstate of at . 

As shown in Refs. lUsl [191 12511 . when the dimension of f) 
satisfies M < poly(A^), the corresponding GCSs play a de- 
cisive role in the theory of entanglement and quantum and 
classical simulations of many-body systems. An example is 
given by the GCSs defined via 

1^/(0) =e-*^^*|0i02---0jv), (23) 

where Hj is the Hamiltonian corresponding to the exactly- 
solvable one-dimensional anisotropic Ising model in a trans- 
verse magnetic field and periodic boundary conditions: 

N 

In sectionliniwe will discuss this system in more detail. 

Any GCS is uniquely determined (up to a global phase) 
by the expectation values of the operators in (j. The state 
|hw(t)) = e~*^*|hw), with if G (), is the highest-weight state 
of () in a rotated CW basis, and satisfies 

/ife(0|hw(0> = Wfc|hw(i)), k G [l,r], (25) 

where hk{t) = e-'"%ke'"' - h + i[hk,H] + ••• G t). 
Thus, 

PhwW - |hw(t))(hw(t)| = J] (/ifc(t) -41), (26) 

where k = W^i^^ivk — u\.) is a constant for normaliza- 
tion purposes. For a particular value of t, the operators 



^k{i) = Sjli ^jit)Qj '^^n be obtained on a classical com- 
puter [i.e., the coefficients _2j{t)] in time polynomial in M 
(see Theorem 1 in Ref. [1^). To see this, note first that 
Xj{t) oc Tr[hk{t)Qj]. Such a trace can be efficiently eval- 
uated by working in the (A/ x M) -dimensional matrix repre- 
sentation (or any other faithful representation) of f) rather than 
working in the (2^ x 2^)-dimensional original representa- 
tion. Therefore, the fidelity of having prepared |hw(i)) can 
be obtained by measuring the expectations of the observables 
appearing in Eq. ( I26I I. over the actual prepared state pi. 

In analogy to the previously discussed cases, a lower bound 
for the fidelity can be obtained by defining the operator 
5ho(t) = [-e(Efe hit) - Vklf] + 1, with e > a con- 
stant determined by the spacing between the eigenvalues ul, 
(see below). If u], < Vk Vz G [1,2^ — 1] one can consider 
S^„{t) = [-£(Efc hkit) - Vkt)] + X instead. Then, 

5h,(t)|0,(t)) (27) 

where |(/!)i(t)) = e^'-^*|(/)i) are the weight states in the rotated 
CW basis (e.g., |hw(t)) = |0o(t))), G K, and = 1. 
Thus, £ is chosen such that Wi = [~s{J2k ■"fe~"fe)^+l] (when 

Vk) + l] (when5„„(t) = [-e(Efe ^fc + satisfies 

Wi<OVi^O. (28) 
For a particular value of t, Eq. ( I28> yields 

^,'„p,„(,) > ('5h.,W)p,. (29) 

This lower bound can be obtained experimentally by measur- 
ing the expectation values of the operators hk{t)hk'{t) and 
(or) hk{t), which are directly induced from the expectations 
{Qj)pi and {QjQj')pi G [l,Af] (assumed to be mea- 

surable with our quantum device). If M = poly(iV) (e.g., an 
evolution due to the Ising Hamiltonian Hj), Eq. j29t can be 
efficiently estimated with C'|poly(Af)] measurements. 

III. QUANTUM SIMULATIONS WITH TWO TRAPPED 
IONS 

In this section we use some of the results obtained in 
Sec. lIICI to estimate the fidelity of evolving two trapped ions 
(qubits) with the Ising-like interaction 

Hi = Jalal + Bial+al), (30) 

where J is the spin-spin coupling and S is a transverse mag- 
netic field. To do so, we will model the system of two 
ions confined in a linear Paul trap and interacting with res- 
onant and non-resonant laser fields as described in Ref. 0]. 
We will estimate the reliability of having prepared the state 
|V>(t)) = |hw(<)) = e-^-^^'lOiOa) (for fixed t). 

In this case, the interaction Hamiltonian for the ions in the 
trap is given by 

Htrap = -ffphonon + -ff|-ionl + -ff|-ion2 + Hm , (31) 
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-ffphonon = '^cmajm'^cm + ^brObr'^br ' 

-ffl-ionl = -[77cmWcm(aL + Ocm) + '7brt^br(abr + abr)]^'^ ' 
H\-\or^2 = -[77cm^cm(aL + Ocm) ^ '7brt^br(abr + "br)]o'a i 

if. = B{al+al). 

Here, the operators (Ocm) ^'I'l "^L ('^br^ create (annihilate) 
an excitacion in the center of mass and breathing modes, re- 
spectively. The coupling interactions i7|_ioni and H\-\ou2 are 
due to the action of state-dependent dipole forces, which are 
generated by the interaction of non-resonant laser beams with 
the electronic levels of the ions (see Ref. 0]). Hm is due to 
the action of an effective magnetic field that can be external 
or generated by resonant laser beams, i/phonon is the energy 
of the normal modes with frequency ti;cm/27r for the center of 
mass mode, and ijjbr/27r for the breathing mode. In the case 
of a single well potential in one dimension, ojbr = \/3i-^cm- 
The couplings (displacements) 77cm and 77br are assumed to be 
small: rj\ <^ 1. They depend on the intensities of the laser 
beams and are given by 

F rir~ 

with i = [cm, br], F the dipole force acting on each ion, and 
m the mass of the ion. 

Therefore, for a fixed value of t, the actual two-qubit state 
prepared in the ion-trap device is 

Pl{t) - Trphonon[e-'^-^*P(ion-phonon)e'^'-^*], (33) 

where we have traced out the vibrational modes. Here, the 
initial state is P(ion-phonon) = |0i02)(0i02| ® Pphonon, and 

Pphonon e T<T^ is the density operator for the initial state 
of the phonons, with the ion motion in a thermal distribution 
being at temperature T (K is the Boltzmann constant). The 
fidelity of having prepared the state \hvj{t)) is then given by 

T\piit),pi,„{t)) ^ Tr[p,(t)phw(t)], (34) 

where the trace is over the spin (i.e., two-qubit) degrees of 
freedom. 

Following the results of Sec. Ill CI we first identify the set 
f)D = Wl,<^z} as the largest set of commuting observables 
in [). This determines |hw) = IO1O2) according to Eq. J21l i. A 
bound for the fidelity of Eq. (|33 can be obtained by using the 
time dependent symmetry operators 

ant) = e-^^'^ale^""'' (j = l,2), (35) 

that uniquely define the state \hw{t)) through the equations 

CT^(t)|hw(t)) = +l|hw(t)). (36) 

Choosing e = 1/2 (see Sec. El} and considering that vi = 
V2 ~ 1, we obtain 51,^^ (t) = ^ [al (t) + al (t)], which satisfies 
[Eq. m] 

^'\pi{t),pUt)) > {Si,nit))pdt)- (37) 



The d-l{t) = {ai — it[Hi,crl] + • • •) are linear combi- 
nations of operators belonging to the Lie algebra 5o(4) = 
{crj, cr^, cr^cr^, o-^CTy , tr^cr^, (y^o-y}- To obtain the coefficients 
involved in these combinations one needs to find the trace be- 
tween the corresponding operators. For example, to obtain the 
coefficient Xi{t) that accompanies the operator al in the de- 
composition of (t), one needs to compute iTr[(T](T](t)]. 
Remarkably, such a trace can be efficiently computed by 
working in the (27V x 2 A^) -dimensional fundamental ma- 
trix representation of so(2A^) rather than in the (2^ x 2^)- 
dimensional original representation (see Ref. fl^ for details). 

In brief, only six correlations (i.e., the elements of so(4)) 
need to be measured to evaluate the inequality of Eq. ( I37l i. 
The complexity of estimating the fidelity is then reduced since 
a naive approach to fidelity estimation would involve the mea- 
surement of fifteen correlations (i.e., the elements of the al- 
gebra su(4)). Of course, the complexity of the problem is 
slightly reduced in this case but the difference is much greater 
for larger systems. 

In Fig. [1] we plot J^^(p;(i), Phw(i)) [Eq- <|34|i] as a func- 
tion of time and for certain values of F, uj\, and B that could 
be attained experimentally. We observe that, for these pa- 
rameters, the fidelity remains close to one, implying that the 
ion-trap device can be used to perform a quantum simulation 
governed by the Ising-like Hamiltonian of Eq. i30\ . We also 
plot {t))pi (t) and we observe that this lower bound of the 
(squared) fidelity already describes much of the reliability of 
the simulation. For the sake of comparison, we also plot the 
expectations (a'J)p^^(t) and {al)p^(ty Finally, inFig.|2]weplot 
the coefficients Xj {t), j £ [1, 6], that determine the weighting 
of the six correlation measurements that contribute to the es- 
timate of (5(,o(i)>p,(t). 



IV. STATISTICAL CONTRIBUTIONS TO MEASURED 
LOWER BOUND ON THE FIDELITY 

In an actual experiment, expectation values can never be 
exactly obtained due to quantum projection noise. Thus, they 
must be estimated after a (typically large) sequence of projec- 
tive measurements performed on identically prepared copies 
of the system. Commonly, maximum-likelihood methods 
(MLMs) I25I27I1 are used to estimate the most probable den- 
sity matrix pi from these measurements. As with full QST, 
these methods are usually inefficient, and they require input 
data concerning every correlation in the system. For example, 
if a MLM is used to estimate the density operator pi of an A^- 
qubit system, the estimation Ep, of the expectation of a par- 
ticular operator E = cr^^ • • • a^^ will require C'[(4^ - 1)A:] 
identically prepared copies of pi, where X is the number of 
copies used to measure a particular correlation (product of 
Pauli operators) j23]- Such a complexity would then be trans- 
lated to the estimation of the lower bounds of Eqs. il2i . jl8l l. 
and ( I29> . In this section we argue that to estimate these lower 
bounds with certain (fixed) level of confidence, the exponen- 
tial complexity can be avoided. 

To prove this, we use results regarding the binomial dis- 
tribution 12^ . Observe first that the operator E, as defined 
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FIG. 1: Numerical simulation of the quantum evolution of two 
trapped ions interacting with laser fields. The parameters used are 
c^cm = WOkHz, 77cm ~ 0.063, B = 560Hz, -J = 5A0Hz, 
F = 25.10"^^iV, and T = 0. These are expected to be attained 
experimentally, (a) Squared fidelity (probability) of having prepared 
the state |hw(t)) = e"'^^* IO1O2), if the dynamics of the trapped 
ions is dominated by the trap Hamiltonian //trap [Eq. <31H . and the 
corresponding lower bound (5i,j-, (t))p,(t), as given by Eq. <37t . as 
a function of time, (b) Expectations of the Pauli operator a], as a 
function of time, if the evolution is governed by Hi and //trap, re- 
spectively. 



FIG. 2: Coefficients \j{t), where alit) = \\{t)Gl + A?(t)cr^ + 
\l{t)(jlol + Xi(t)olGl + Af (f)(J>2 + A?(t)(j>2 £ so(4)^ used 
to obtain (<5i,o)p, (t) in Fig.Q Note that, because of the symmetry 
under permutation of both ions, the same coefficients are obtained in 
the decomposition of cr^{t). 



Eq. ( I38t guarantees that Sp, differs by at most \/TJX from 
the actual expectation with (at least) 68% of confidence I30I1 . 
For example, if S is estimated from ten thousand identical 
copies of pi, then (E)p, = Sp, ± .01 with (at least) 68% of 
confidence. 

With no loss of generality, the bounds of Eqs. \\2\ . \\%\ . 
and ( I29> . can be rewritten as 



^ > ao + 



R 
m— 1 



j(I]"^)pj , ao, am € 



(40) 



where each S™ involves a particular product of Pauli oper- 
ators {R ~ poly(A^)]. If each is estimated from X 
identical copies of pi, then (E™) 
of confidence, and 



/Pi 



± V'T7^with68% 



above, has ±1 as possible eigenvalues. Then, if we perform 
projective measurements of S over X identical copies of pi, 
we obtain 

(S)p, = Ep, ± 5, (38) 

where Ep^ = is the estimated expectation (i.e., X± 

are the number of times we measured E = ±1, respectively), 
and 5 is the corresponding standard deviation. The latter is 
given by 

5 = 2^^, (39) 

where p± are the (not known) probabilities of measuring E = 
±1, respectively. Then, 6 < yfXjX. 

For sufficiently large X, the binomial distribution can be 
well approximated by the normal distribution. In this context. 



li R 
rn— 1 m— 1 

(41) 

with the same confidence. Of course, Eq. ( I41> provides rel- 
evant information if E™ 3> \j\fX. For example, if one 
is interested in preparing the state IGHZtv), then R = N 
and Ep" w +1. Choosing X = 10'' iV'^, a good estimation 
(with error 0.01) for the lower bound of the fidelity is ob- 
tained. The method is then efficient: lower bounds on fidelity 
of state preparation can be obtained, with certain confidence, 
in poIy(iV) identical preparations of pi. 

We have not considered any source of error other than the 
one given by the statistics of projective measurements in the 
quantum world. Otherwise, the results obtained in the pre- 
vious sections must be modified according to the specific 
sources of error or decoherence that can affect the state prepa- 
ration. 
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V. CONCLUSIONS 

We have studied the fidelity of state preparation for three 
different classes of states: the rotational-invariant states, SSs, 
and GCSs. Many interesting multi-partite entangled states, 
like cat or W-type states, belong to these classes. In partic- 
ular, GCSs are natural in the framework of quantum simula- 
tions. We have discussed the quantum simulation of the two- 
qubit Ising model using an ion-trap device. In this case we ob- 
serve that a lower bound of the fidelity of the simulation can 
be simply obtained, and can be considered to estimate the reli- 
ability of the experiment. Such a bound can also be efficiently 
estimated for other multiple qubit systems having Ising-like 
interactions. Similar approaches can be considered to study 
the fidelity of state preparation in general qudit or fermionic 
systems. 

Our results provide an efficient method to estimate, with 
certain confidence, lower bounds on the fidelity of state prepa- 



ration based on symmetries. Many of the states described 
contain A^-particle entanglement, so the lower bounds can 
also be used to verify entanglement using entanglement wit- 
nesses I31I l32il . These bounds are most accurate when the 
actual prepared state is not too far from the desired one, as in 
Fig E Therefore, a consequence of our results is that instead 
of measuring every possible quantum correlation of a system 
a large number of times (as for QST), one should focus on 
having good estimations of certain relevant expectations. 
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